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Charmonium production from nonequilibrium charm and anticharm quarks in 

quark-gluon plasma 

Taesoo Song/'E Kyong Choi Han, 2 -0 and Che Ming Ko 2 @ 

'Cyclotron Institute, Texas AhM University, College Station, Texas 77843-3366, USA 
2 Cyclotron Institute and Department of Physics and Astronomy, 
Texas Ak,M University, College Station, Texas 77843-3366, USA 

Parameterizing the charm and anticharm quark momentum distributions by the Tsallis distri- 
bution, we study the nonequilibrium effect on the charmonium production rate in a quark-gluon 
plasma up to the next-to-leading order in perturbative QCD. We find that nonequilibrium charm 
and anticharm quarks suppress the charmonium production rate compared to that from equilibrated 
ones. We further show that the suppression factor calculated with the charm quark relaxation time, 
which has been frequently used in the literature, is close to our results. 

PACS numbers: 25:75.Cj 



I. INTRODUCTION 



Since the suggestion of using charmonium suppression 
as a possible signature for the formation of a quark-gluon 
plasma (QGP) in relativistic heavy-ion collisions [l|, 
there have been many theoretical |7H15j and experimen- 
tal [2H6( studies on charmonium production in heavy ion 
collisions at energies available from the Super Proton 
Synchrotron (SPS), the Relativistic Heavy Ion Collider 
(RHIC), and the Large Hadron Collider (LHC). Although 
the collision energy spans two orders of magnitude from 
SPS to LHC, the observed suppression of charmonium 
production in these collisions was found to be similar. 
A possible explanation for this surprising result is that 
the expected larger suppression of charmonium produc- 
tion with increasing collision energy is compensated by 
the increasing number of regenerated charmonia from ini- 
tially produced charm and anticharm quarks in QGP. 
This is because the number of reg enerated charmonia is 
quadratically proportional to [16J ] while that of directly 
produced charmonia depends linearly on the number of 
charm quarks J16H18J . Studies based on phcnomenologi- 
cal models indeed show that although the fraction of re- 
generated charmonium is compatible to that of directly 
produced one in central Au+Au collisions at RHIC, it be- 
comes dominant in central Pb+Pb collisions at LHC [lOl — 
Il4| as more charm and anticharm quarks are produced. 

The charmonia production rate in QGP depends on 
the charm and anticharm quark distributions in both 
momentum and coordinate spaces [l9j . Presently it is 
not clear to what degree charm and anticharm quarks, 
which arc produced from initial hard collisions with a 
power law distribution, are thcrmalizcd in relativistic 
heavy-ion collisions. It seems that charm and anticharm 
quarks of low transverse momentum are close to while 
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those of intermediate and high transverse momenta are 
far from thermal equilibrium [2CJ|. In Refs. [Hi - Ha l . the 
charm quark nonequilibrium effect on charmonia regen- 
eration in QGP is included by multiplying the expected 
equilibrium charmonium number by a relaxation factor 
R(t) = 1 — cxp[— (r — To)/r re i] with r re i being the re- 
laxation time of charm quarks in QGP. Since the latter 
has a value T rc i ~ 4 fm/c, the regeneration contribution 
to charmonium production is suppressed by a factor of 
0.53, 0.68 and 0.8 in heavy ion collisions at SPS, RHIC, 
and LHC, respectively, when the corresponding QGP life- 
time is about 3.0, 4.5, and 6.5 fm/c. To check the ac- 
curacy of the estimate based on the relaxation factor, 
we calculate in the present study the charmonium pro- 
duction rate using the Tsallis distribution for charm and 
anticharm quarks [2l|, 122] ■ The Tsallis distribution with 
a parameter q is a generalized Boltzmann distribution, 
which is the same as the normal Boltzmann distribution 
when q = 1 but becomes a harder power law distribu- 
tion as the value of q increases. In this study, we define 
A = q — 1 and determine its values for the initially pro- 
duced charm and anticharm quarks from hard collisions 
at RHIC and LHC by fitting the charm quark spectra 
generated by the event generator PHYTIA ,23j . To study 
how the charm and anticharm quarks approach equilib- 
rium in a QGP, we use the parton cascade model based 
on their elastic scattering with light quarks and gluons in 
the QGP. Fitting these distributions by the Tsallis dis- 
tribution, we then calculate the charmonium production 
rate using the transition amplitude that is calculated up 
to the next-to-leading order (NLO) in perturbative QCD 
(pQCD) [24|,[25||. We find that nonequilibrium charm and 
anticharm quarks suppress the charmonium production 
rate compared to that from equilibrated ones and that 
the suppression factor estimated from the charm quark 
relaxation time, which has been frequently used in the 
literature, is close to our results. 

This paper is organized as follows. We first give in 
Sec.[TT]the transition amplitudes for charmonium produc- 
tion in pQCD. In Scc. lIIIl we derive the expression for the 
charmonium production rate from charm and anticharm 



quarks based on the Tsallis distribution and study the 
dependence of the J/tp production rate on the parame- 
ters in the Tsallis distribution. We then apply in Sec. IIVI 
these results to heavy-ion collisions at RHIC and LHC. 
Finally, we give the conclusions in Sec. |Vj For details on 
the evaluation of the phase space integral in the charino- 
nium production rate, they are given in Appendix A. 



II. TRANSITION AMPLITUDES FOR 
CHARMONIUM PRODUCTION 

Although the suggestion of J/tp suppression as a sig- 
nature of QGP was based on the idea that the J/tp could 
not be formed in QGP [1| due to screening of the color 
charge, recent studies on the spectral functions of char- 
monia have suggested, on the other hand, that the dis- 
sociation temperature of J/tp is higher than the critical 
temperature for QGP phase transition |26l - l28j . indicating 
that the J/tp can survive and be regenerated in QGP. 

In pQCD, the leading-order process for J/tp production 
from charm and anticharm quarks is the reaction c + c — > 
J/tp + g. The transition amplitude for this reaction is 
same as that for the J/t p dissociation reaction J/t/j + g — ► 
c + c given in Rcfs. [2J, [25J, that is 
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In the above, fco and p = \Jm c (k$ — eo) are, respectively, 
the gluon energy and the relative three momentum of 
charm and anticharm quarks in the J/tp rest frame with 
eo being the binding energy of J/tp; g is the strong cou- 
pling constant which is taken to be 1.87 from previous 
phenomcnological studies (14J |29|. and m c and mj/^ are 
the masses of charm quark and J/tp, respectively; m g is 
the thermal gluon mass obtained from lattice data using 
the quasiparticle model [30(; and tp(p) is the wavefunc- 
tion of J/tp. To include the medium effect, we use the 
temperature-dependent wavefunction and binging energy 
of J/tp obtained from solving the Schrodinger equation 
using the screened Cornell potential between charm and 
anticharm quarks |3l|. Eq. ([I I) becomes the same as that 
of Bhanot and Peskin [32|, [33| if tp(p) is taken to be the 
Coulomb wavefunction for the 15* state. 

For the next-leading order (NLO) reaction c + c + 
Q(Qi9) ~^ J/'ip + <z(<Zj<?) m pQCD, its transition ampli- 
tude is the same as that for the reaction J/tp + q(q, g) — > 
c + c + q(q,g). As given in Ref. J2J, [25(, the squared 
transition amplitudes of quark- and gluon-induccd NLO 
processes are 
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where N c is the number of quark colors, and fci and k 2 
are the momenta of incoming and outgoing partons, re- 
spectively. 

Other NLO reactions include c + c^ J/4> + g + g(q + q) 
which can be obtained from the reactions c+c+q{q, g) — > 
J/tp + q(q, g) by changing the incoming parton to an out- 
going parton. Their squared amplitudes can thus be ob- 
tained from Eqs. ^ and ^ by changing fci and fcio to 
—fci and — fcio, respectively, and by multiplying Eq. ^ 
by an overall minus sign and the nubmer of light quark 
flavors, that is 
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III. THE CHARMONIUM PRODUCTION RATE 

The production rate of charmonium with momentum 
q in a QGP from the LO and NLO charm and anticharm 
quark recombination processes described in the previous 
Section can be written as u3 
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In the above, / c (pi) and / e (P2) are, respectively, the 
charm and anticharm quark distributions; / p (fci) is the 



thermal parton distribution including the degeneracies 
of quarks due to their flavors and antiparticlcs; k is the 
outgoing gluon momentum in the LO c + c -+ J/tp + g 
reaction; and k.2 is the outgoing parton momentum, while 
ki is the outgoing and incoming parton momentum in the 
NLO c + c ^ J/tp + g + g(q + q) and c + c + q(q, g) -> 
J/tp + q{q, g) reactions, respectively. 
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FIG. 1: (Color online) Charm quark distribution for different 
values of A in the Tsallis distribution. 
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FIG. 2: (Color online) Production rate for charmonium of 
momentum q in QGP at temperatures T = 1.0 (left panel) 
and 1.5 T c (right panel) from charm and anticharm quarks 
that have a Tsallis distribution with various values of A. 



charm quark momentum spectrum. Also, the J/tp pro- 
duction rate is larger at 1.5 T c than at 1.0 T c , consistent 
with the larger J/tp thermal decay width at higher tem- 
perature. 



For the charm and anticharm quark distributions, we 
use the modified Tsallis distribution function 
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where A = q — 1 with q being a parameter in the original 
Tsallis distribution function and A, which depends on A, 
is given by 
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to ensure that the integration of / c (p) over the charm 
quark momentum gives the charm quark density n c . In 
Fig.Q] we show the charm quark distribution for different 
values of A. It is seen that as A approaches 0, f c becomes 
the Boltzmann distribution. Details on the phase-space 
integrations in Eqs.©, (JT]), and © with the Tsallis dis- 
tribution are given in Appendix A. 

In Fig. [2j we show the production rate for charmo- 
nium of momentum q in QGP at temperature T = 1.0 
(left panel) and 1.5 T c (right panel), where T c is taken to 
be 170 MeV, from charm and anticharm quarks that have 
a Tsallis distribution with various values of A =0, 0.04, 
0.08, and 0.12. These results are obtained with the den- 
sity of charm quarks in the QGP taken to be 0.008/fm 3 
(1 GeV 3 ). It is seen that the momentum dependence of 
the J/tp production rate becomes harder as A increases. 
This is reasonable because a larger A indicates a harder 
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FIG. 3: (Color online) Production rate of J/tp per unit rapid- 
ity as a function of A at T — 1.0 and 1.5 T c . Dotted, dashed, 
and dot-dashed lines are from the reactions c + c — s> J/tp + g, 
c + c + q(q,g) -»■ J/ip + q(q,g), and c + c -+ J/tp + g + g(q + q), 
respectively, and solid lines are their sum. 

Figure [3] shows the dependence of the momentum in- 
tegrated J/tp production rate per unit rapidity on the 
parameter A for the three reactions c + c — > J/tp + g (dot- 
ted line), c+c+q(q, g) -+ J/tp + q(q,g) (dashed line), and 
c+ c — >• J/tp + g + g(q + q) (dot-dashed line). The most 
dominant process is c + c + q{q, g) — > J/ip + q(q, g), es- 



pccially at high temperature. The reason for this is that 
the common factor |9-0(p)/9p| 2 in the squared transi- 
tion amplitudes in Eq. ([I])—© peaks at small relative 
momentum p between charm and anticharm quarks and 
the value of p is smaller in the reaction c + c + q(q, g) — > 
J/ip + q{q, g) than in other reactions. 

Figure [3] also shows that the J/ip production rate from 
the three reactions behave differently as the value of A 
increases. Since harder charm and anticharm quark spec- 
tra are less favorable to J/ip production as a result of the 
large relative momentum between charm and anticharm 
quarks, the production rate of J/ip from the reaction 
c+c+q(q, g) — > J/ip+q(q, g) decreases as A increases, sim- 
ilar to that found in the coalescence model QJl . Although 
the transition amplitudes in all processes favor soft charm 
and anticharm quark spectra, to produce additional par- 
tons in the final state as in the reactions c + c^ J/ip + g 
and c + c ^ J/ip + g + g(q + q) requires that the initial 
charm and anticharm quarks to have high energy. As a 
result, the production rates of J/ip from the reactions 
c + c^)- J/ip + g and c + c^>- J/ip + g + g(q + q) increase 
and then decrease with increasing A, particularly from 
the latter reaction as it has two thermal partons in the 
final state. The total J/tp production rate given by the 
sum of the production rates from the three reactions de- 
creases, however, with increasing A since it is dominated 
by the reaction c + c+q(q, g) — >■ J/ip + q(q,g). Therefore, 
noncquilibrium charm and anticharm quark distributions 
suppress the J/ip production rate in QGP compared to 
that from completely thermalized charm and anticharm 
quark distributions. 



IV. APPLICATIONS TO HEAVY-ION 
COLLISIONS 

To see the relevance of our results to J/ip production 
in heavy ion collisions, we carry out a cascade calculation 
to study how the charm and anticharm quark distribu- 
tions change with time in a QGP at a fixed temperature. 
Specifically, we take the initial charm and anticharm 
quark spectra at mid-rapidity (\y\ < 1) from the PHY- 
TIA simulations [23[ for p+p collisions at yfs = 200 GeV 
and 2.76 TeV, available at RHIC and LHC, respec- 
tively, and they are shown by solid and open circles in 
panel (a) of Fig. 0J Assuming that the charm and an- 
ticharm quark spectra are isotropic in momentum space, 
we then find that they can be very well described by the 
Tsallis distribution with (T, A) =(650 MeV, 0.013) and 
(680 MeV, 0.033) as shown, respectively, by lines in panel 
(a) of Fig. [4] Taking the elastic scattering cross sections 
of charm and anticharm quarks with light quarks and an- 
tiquarks to be 1 mb and those with gluons to be 2 mb, 
we examine the change of their momentum distributions 
with time. The results for the two initial charm and an- 
ticharm quark spectra corresponding to \/s ~ 200 GeV 
and 2.76 TeV in a QGP of temperature T = 1.5 T c are 
shown by solid and open circles in panels (b)-(d) of Fig.0] 
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FIG. 4: (Color online) Charm quark momentum distribu- 
tion in arbitrary units for RHIC (solid circles) and for LHC 
(open circles) in QGP of T = 1.5 T c at 0, 1, 2, and 3 
times the charm quark relaxation time t Te i. Lines are fit- 
ted Tsallis distributions with (T, A) =(650 MeV, 0.013), 
(190 MeV, 0.078), (140 MeV, 0.072), and (150 MeV, 0.058) 
for RHIC and (T, A) =(680 MeV, 0.033), (140 MeV, 0.107), 
(90 MeV, 0.097), and (100 MeV, 0.081) for LHC at increasing 
time. 



for different times of li ro i, 2t ro i, and 3i re i. Here t re \ de- 
notes the charm (anticharm) quark relaxation time and 
is defined by 
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where rii is the density of parton species i includ- 
ing its degeneracy factor in grand canonical ensemble, 
v Te i is the relative velocity between the charm (an- 
ticharm) quark and parton, Ui is the elastic cross sec- 
tion of charm (anticharm) quark by the parton, p and 
p' are three momenta of incoming and outgoing charm 
quarks in laboratory frame, respectively, and (• • •) de- 
notes an average over the charm quark momentum dis- 
tribution. The relaxation time is thus the inverse of 
drag coefficient of charm quarks [35j. We find from 
the parton cascade simulation that the relaxation time 
of charm quarks in a QGP of temperature 1.5 T c de- 
pends only weakly on their momentum distribution and 
has an approximate value of about 3.2 fm/c. In Fig. 0] 
we also show by lines the fitted charm (anticharm) 
quark momentum distributions using the Tsallis distri- 
bution with the parameters (T, A) = (190 MeV, 0.078), 
(140 MeV, 0.072), and (150 MeV, 0.058) for the case 
of RHIC and (T, A)=(140 MeV, 0.107), (90 MeV, 0.097), 
and (100 MeV, 0.081) for the case of LHC at different 
times of l£ re i, 2i re i, and 3i rc i. It is seen that the Tsal- 
lis distribution again describes very well the charm (an- 
ticharm) quark distribution from the cascade simulation. 
We note that the results shown in Fig. @] with the time 



measured in unit of the relaxation time hardly depends 
on the value of elastic cross sections of charm and an- 
ticharm quarks. 



quarks, the J /if) production rates in heavy ion collisions 
at RHIC and LHC are then suppressed, respectively, by 
40-50 % and 35-45 % compared to those from completely 
thcrmalizcd charm and anticharm quarks. 
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FIG. 5: (Color online) Time dependence of the ratio of the 
J/ip production rate from nonequilibrium charm and an- 
ticharm quarks to that from completely thermalized charm 
and anticharm quarks in QGP at T = 1.0 and 1.5 T c . The 
solid line is the result based on the relaxation factor correc- 
tion. 

Using the J/ip production rate given in Fig. [3] and the 
charm and anticharm quark distribution given in Fig. |4l 
we can calculate the J/ip production rate and obtain its 
ratio with respect to the rate from completely thermal- 
ized charm and anticharm quarks. This is shown in Fig.[S] 
as a function of time for a QGP of temperature T =1.0 
and 1.5 T c in collisions at RHIC and LHC. These re- 
sults show that the suppression of J/ip production due 
to the nonequilibrium charm and anticharm quark dis- 
tributions is stronger at higher-energy collisions and at 
higher temperature. This is due to the fact that the 
more energetic charm and anticharm quarks produced in 
higher energy collisions make it harder to produce the 
charmonium, and that the smaller charmonium binding 
energy at higher temperature further suppresses its pro- 
duction. Although the suppression starts at different val- 
ues in different cases, their time dependences all seem 
to be similar. Since the temperature of QGP formed 
in relativistic heavy ion collisions decreases with time, 
the suppression of J/ip production due to nonequilib- 
rium charm anticharm quarks would decrease with time 
as well. For comparison, we also show in Fig. [5] by the 
solid line the results from using the relaxation factor 
R = 1 — exp[— r/r re i] that has been used in the liter- 
ature to take into account the effect of nonequilibrium 
charm and anticharm quark distributions [IJ, [jjl, LL4| . It 
is seen that the suppression of the J/ip production rate 
based on the nonequilibrium Tsallis distribution is sim- 
ilar to the relaxation factor correction. Assuming the 
lifetime of QGP is given by the relaxation time of charm 



V. CONCLUSIONS 

Using the transition amplitudes for J/ip production 
from charm and anticharm quarks calculated up to NLO 
in pQCD and including the medium effect on the J/ip 
wavefunctions and binding energies from a screened Cor- 
nell potential between the charm and anticharm quarks in 
QGP, we have studied charmonium production in QGP 
from nonequilibrium charm and anticharm quarks that 
are described by the Tallis distribution. We have found 
that nonequilibrium charm and anticharm quarks sup- 
press the production rate of J/ip in QGP, compared 
to the rate from completely thermalized charm and an- 
ticharm quarks. We have further used the calculated 
J/ip production rate in QGP to study J/ip production 
in heavy-ion collisions at RHIC and LHC by using the 
PHYTIA simulation to obtain the initial charm and an- 
ticharm quark distributions and then follow their clastic 
scattering with light quarks and antiquarks as well as 
gluons in QGP via the cascade simulation. With the re- 
sulting charm and anticharm quark distribution param- 
eterized by the Tsallis distribution, we have found that 
the suppression in the J/ip production rate is stronger in 
higher energy collisions and at higher temperatures, al- 
though their time dependence is similar in all cases. We 
have also compared our results to those in the literature 
obtained with the relaxation factor correction using the 
charm quark relaxation time and found that the latter is 
similar to our results. 
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Appendix A 

In this Appendix, we give details on the evaluation 
of the J/ip production rates in QGP from the reactions 
c + c^ J/ip +g, c + c^i- J/ip + g + g(q + q), and c + c + 
Q.i.0.1 9) ~^ J/ 1 !' + Q.i.0.1 9)1 based on the Tsallis distribution 
for the charm and anticharm quark distributions. 



c+ c 
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C + C + q(q, g) -> J/tp + q(q, g) 



Substituting Eq. (|9|) into Eq. (]6j) and using d 3 q = 
EqirnTdrnxdy in terms of the J/ip transverse mass rrvr 
and rapidity y, we obtain 
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where V is the volume of the system. To evaluate above 
integrals, we first express the Tsallis distribution function 
in covariant form by replacing E p with p ■ u, where u M = 
7(1, 0) with p being the boost velocity and 7 = (1 — 
/3 ) ' . Transforming to the center of mass frame of 
charm and anticharm quarks via P = Pi + P2 and p = 
(pi - Pa)/2, which leads to j3 = V jE and 7 = E/y/s, 



with E = E p 

E 2 — P 2 , we then have 



E P2 and s 



(q + k) 2 = [ P i+p 2 f 



As in the previous case, the production rate from this 
reaction can be expressed in the center of mass frame of 
charm and anticharm quarks, and it is given by 



diVNLOl 



A 2 (\)V 
rriTdrriTdydt 16(2-7r) 3 



d 3 ki 



d 3 k 2 



1-^\M\ 

w 



2 
NLO 



(2n)32E kl J (27r)32£ k2 
1 
i/( k i) / 



-1 
AP\ 2 



A7v / w\ 2 /, Am 2 * 

-&-) -\ 1 -—) 1 2T ; cos " 



/AFy 



-l/A 



(A4) 



where ki and k2 are, respectively, the three momenta 
of incoming and outgoing partons; P = pi + P2 and 
w = (pi + p 2 ) 2 ; and /(ki) is the Fermi-Dirac distribu- 
tion in the quark-induced reaction and the Bose-Einstein 
distribution in the gluon-induced reaction. Furthermore, 
/(ki) should be multiplied by the flavor number in the 
case of quark-induced reaction. 

Choosing the z-axis along the J/tp momentum and us- 
ing 



dN LO 



mTdrriTdydt 



A 2 (X)V 
8tt 2 



d 3 k 



d 3 p 

(2n) 3 2E w J (2n) 2 s 



(} + *£)'-( 



^(Vs-E^-E^Mlto 



A7/3-p\ 2 



T 



-l/A 



d 3 k 



^ 2 (A)y 

16(2tt) 3 J (2n) 3 2E k 



Am 2 



\M\l 



dcoa9 



Ami 



AP\2 
V27V 



(— Y C os 2 e 



-l/A 



(A2) 



where is the angle between p and (3. 

Choosing the momentum q of produced J/tp to be 
the z-axis, i.e., q^ = (-E q ,0,0, g), and denoting k v = 
(-Ek, 0, k sin </>, k cos (p), we then have 



dNff* 



A 2 (X)V 
mTdrriTdydt 32(27r) 5 



dkk 2 
~Ek~ 



1-M| M | 



2 

LO 



dcosi 



dcos9 



1+^) -U- 



AP 
2T 



-l/A 



k" = (-Eki, 0, k\ sin0, k\ cos</>), 
k 2 = (E^, k 2 sin 9i cos 9 2 ,k 2 sin 9i sin 9 2 ,k 2 cos 9i), 
P° = {E^ + E^-E^, -k 2 sm6 1 cos0 2 , 

k\ sin (p — k 2 sin 9i sin 2 , k\ cos <p — k 2 cos 9\ ) , ( A5) 

where P = p\ + p 2 , the production rate becomes 
A 2 (\)V f°° dkxkl rb 



jAfNLOl 



rriTdrriTdydt 64(27r 
^2 



dk 2 k? 



E k2 



E kl 
d cos 0i 



dcos( 



i-M|^| 



2 

NLOl 
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where 
P 2 

s 



A7\/w \ 2 / 4 ? 7t 2 
2T ) \ w 



1 
d0 2 

dcos9 



/AP\2 ' 



-l/A 



(A6) 



= q 2 + k 2 +kj + 2g(fciCOS(/)- fc 2 cos0i) 
— 2fcifc 2 (cos</i>cos0i + sin0sin0i sin0 2 ), 



(Eq + -Eki ~~ ^k 2 



V/V-' 



l ?(?,s) 



2E< l E kl - 2gfci cos 0. (A7) 



, Because the incoming energy is always larger than the 

invariant mass of final state, i.e., s > (2m c + "i g (g, g )) 2 ; 
(A3) fei and cos have the following integration limits 



where P 2 = q 2 +k 2 + 2qk cos <p. We note that the squared 
transition amplitude, Eq. ([T]), is a function of ko, which 
is the energy of thermal gluon in the J/tp rest frame and 
can be expressed as fco = (s — mL — m 2 )/2mj/ t p. 



-C + \/C 2 - ABD 



v/i/> 



2B 

- Am 2 + Am c m a ( g 



q(q,9) 



2-Eq-Eki 



2qki 



, (A8) 



with 



the J/if> production rate is then 



B = 4m^, 

C = Aq{Am 2 c - m 2 ^ + 4m c m g ( ?iS )), 

D = AElm 2 q{qg) - {Am 2 ~ m j/^ + 4m c m gfeg) ) 2 . 

For the integration ranges of k 2 , cos#i, and 9 2 , they are 
such that w > Am 2 .. 

The squared transition amplitude for this reaction, i.e., 
Eq. ((2J or ([3]), is a function of the energies feio and £20 
of incoming and outgoing partons given by 



ho 

with 



» - m ?W - m s 



2m 



J/V 



«20 



'J/4, 



2m 



J/V 



w = (q- k 2 y 



l q(q,g) 



-2EqE k2 + 2qk 2 cos9 1 , 



(A10) 



dN j/f 2 _A 2 (X)V /-"dkik? 



,2 /•! 



rnxdrnxdydt 64(27r) 8 Jq i? kl j_] 

<.27T 

d cos Q\ I d9 2 



dcos( 



11 



-Ek 2 



-1 



II 



XA 1 



M| M | 



2 
NL02 



dcos^ 



-, (A9) where 



(i+^y-o-^) 



m?\/AP\2 2 
1 — I cos" 2 



\2T) 



-l/A 
(A14) 



q 2 + k 2 + k 2 . + 2g(fci cos <j> + k 2 cos #1 ) 
+2fcifc2(cos</i>cos#i + sin</>sin#i sin 6*2), 
(£ q + £ kl +£ k2 ) 2 -P 2 . 



and of 



h • h = Ski-Eka 

— fcifc2(cos<?!)cos6'i 



smfflsint'i sm tl 2 ) 



(All) 



The squared transition amplitude, i.e., Eq. Q or ([5]), is 
a function of the energies fcio and k 2 o of thermal partons 
in the J/tp rest frame, 



3. c + c^ J/ip + g + g(q + q) 

The production rate for this reaction is similar to the 
one for the reaction c + c ^ J/ip + g, i.e., 



dN*%» 



A 2 {X)V 



d 3 ki 



d 3 k. 



m T dm T dydt 16(2tt) 3 J (27r) 3 2£ kl y (27r) 3 2£ k2 

r-l 

<icos0 



-^|M| 



2 
NL02 



('+¥)'-(' 



4m^ 



\2TJ 



-l/A 

(A12) 



CIO 



where 



s i - m j/0 - TO s 



2m 



J/V 



s 2 - m 2 7/ . - m 2 
k 20 = ^ %M5) 



2m 



J/1> 



8! = (q + hf 

= (£q + E kl f - q 2 - k 2 - 2qk 1 cos < 
s 2 = {q + k 2 ) 2 



= (E q 



E u 



— q 2 — kl - 2qk 2 cos ( 



where h and k 2 are momenta of two outgoing thermal 
partons. 

Again, choosing the z-axis along the J/tJj momentum 
and using 



k\ = (£^ki ! 0, ki sin <f>, ki cos <j>) , 
k£ = (-E k2 , k 2 sin 9 1 cos 2 1 k 2 sin 9\ sin 02 , k 2 cos 9\ ) , 



(A13) 



and of 



h ■ k 2 



S 2 + m j/^ 



(A16) 
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